Unconventional Hall effect in pnictides from interband interactions 
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We calculate the Hall transport in a multiband systems with a dominant interband interaction 
between carriers having electron and hole character. We show that this situation gives rise to an 
unconventional scenario, beyond the Boltzmann theory, where the quasiparticle currents dressed by 
vertex corrections acquire the character of the majority carriers. This leads to a larger (positive 
or negative) Hall coefficient than what expected on the basis of the carrier balance, with a marked 
temperature dependence. Our results explain the puzzling measurements in pnictides and they 
provide a more general framework for transport properties in multiband materials. 
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The discovery of superconductivity in iron-based su- 
perconductors has triggered a renewed interest in the 
properties of interacting multiband systems. Indeed, all 
the families of pnictides display several hole (h) and 
electron (e) pockets at the Fermi level, as predicted 
by density-functional-theory calculations and confirmed 
by Fermi-surface sensitive experiments, as de Haas van 
Alphen and photoemission spectroscopy [1]. A much 
more indirect probe of such a multiband character comes 
from transport experiments, where the contribution of 
carriers having h and e character is unavoidably mixed. 
A typical example is provided by Hall-effect measure- 
ments, for which a standard Boltzmann-like multicarrier 
picture [2] would give a Hall coefficient 
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where n a (a =e, h) is the density of the e- or h- car- 



rier type, and \x a 



~T a /m a the corresponding mo- 



bility, with T a being the transport scattering time and 
m a the effective carrier mass in each band. In compen- 
sated semimetals as pnictides, where n e ~ n^, one would 
thus expect an almost vanishing Hall coefficient. Strik- 
ing enough, most pnictides [3-9] show a quite different 
scenario with a very large absolute value of the Hall co- 
efficient Rh, and with a marked e/h character of the 
transport in materials which are only slightly e/h doped. 
In addition, a strong temperature dependence of Rh is 
also typically found, that disappears only at very large 
doping away from half- filling. Since \i a oc r a , to account 
for these features within the Boltzmann-like approach (1) 
one needs thus to assume a marked disparity between r e 
and Tfi , with r e 3> th in e-doped compounds and th 3> r e 
in the h-doped ones. However, such a disparity has not 
been supported until now by any explicit calculation. For 
example, the inclusion of spin [10] or orbital[ll] fluctua- 
tions within realistic models can account at most for a 
factor of 2 of anisotropy between the average quasiparti- 
cle lifetime on the e and h pockets, not enough to explain 



neither the absolute value of Rh nor its T dependence 
reported in Refs. [3-9] . In addition, the claim of two very 
different scattering rates in pnictides obtained from op- 
tical probes, where the fiat mid-infrared optical conduc- 
tivity is sometimes attributed to a very broad Drude-like 
intraband contribution [12-14], has been questioned by 
several authors on the basis of the presence in pnictides 
of low-energy interband optical transitions [13, 15-17]. 
A convincing framework to explain the unconventional 
properties of the Hall transport in pnictides, and hence 
to elucidate the role of the scattering mechanisms, is thus 
still lacking in these materials. 

In this Letter we show that in a multiband system with 
predominant interband interactions between carriers hav- 
ing opposite (e/h) character the semiclassical picture of 
transport based on Eq. (1) must be strongly revised. By 
computing explicitly the current vertex corrections due to 
the exchange of spin fluctuations between h and e bands, 
we show that, in contrast to the standard Fermi-liquid 
case, they cannot be simply recast in a renormalization 
of the transport scattering time with respect to the quasi- 
particle lifetime. Indeed, the spin fluctuations induce a 
mixing of the electron and hole currents such that the 
renormalized current in each band can even have oppo- 
site direction with respect to the bare band velocity. This 
mechanism explains the large value of \Rh\ in slightly 
e/h doped compounds, and its temperature and doping 
dependence, in good agreement with the experimental 
findings in the non-magnetic state. 

Let us introduce the minimal model which contains the 
main ingredients responsible for the unconventional Hall 
transport in pnictides. We consider a two-band model 
with two-dimensional parabolic e/h bands centered at 
the r and M = Q = (tt, 7t) points, with different Fermi- 
surface (FS) areas: 
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where k is the reduced momentum with respect to the Y 



2 



or M point, for the h or e band, respectively. We take 
here units ft = c = a = 1 (a being the lattice spacing). In 
the following we choose fi — 0, so that E^ in and E^ ax 
fix the Fermi wavevectors k F in each band, and we will 
assume, without loss of generality, m e — mh = m. The 
general gauge-invariant expression for the longitudinal 
and transverse conductivities for each band a:=e,h can 
be derived on the basis of the Kubo formula for a weakly 
interacting system [18, 19]. In particular we have that 
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where the derivative of the Fermi function f(x) = (1 + 
e x / T )~ x has been approximated at low temperature T 
with a 5-function, so that only quantities at the FS ap- 
pear in the final expression. In Eq. (3) the vector v Q 
denotes the band velocity in the x-y plane for the band 
a, 3 a is the corresponding renormalized current and T a is 
the inverse quasiparticle lifetime, determined in general 
by electron-electron and impurity scattering processes. 
Due to the symmetry of the problem, J Q is parallel to 
the reduced moment k in each band a, and r a (k) and 
J"(k) depend only on |k|, so we define their value at the 
FS as T F = T a (kp) and 3 F = 3 a (k F ). The transverse 
xy conductivity under a weak magnetic field H along the 
z axis can be written as: 
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Here e z is the 
= (e, x v a )/|v Q | 



where A a (k) = v a [j Q x (e| • V)J Q 

unit vector along the z axis, while e' 
is tangential to the a-th FS at k. For a parabolic band 
we have thus e" ■ V = ±de/k, where the plus/minus 
sign holds for an e/h band, respectively, so that A a = 
±{v a /k)(3 a x d 3 a ) z = ±v a (J a ) 2 /k. As a consequence, 
the overall sign in Eq. (4) is determined only by the 
sign of v • k, which identifies the e/h character of the 
band [20]. Once evaluated the longitudinal and transverse 
conductivity, the Hall coefficient Rh is given by 



Rh = 



xy 



(5) 



Eqs. (3)- (5) are quite general, since they express the 
conductivities in terms of the bare Fermi velocity v° and 
the renormalized current J". What makes multiband 
systems peculiar is the nature of vertex corrections that 
determine the relation between v Q and 3 a . Using a stan- 
dard approach[21] one can establish between these two 
quantities a matricial relation 



Jf — A a/ 3V F . 



(6) 



The case A a p — 5 a fj corresponds to the non-interacting 
system where, by using the 2D relation n a = (k F ) 2 /2tt 



and by identifying l/r a = 2T F , Eqs. (3)- (4) reduce to the 
standard results, a" x = e 2 n a T a /m and a" y = ^fa^ x /j, a H, 
with [i a = eT a /m, and the minus/plus sign holds for 
the e/h band, respectively. In the interacting case, the 
strength of the diagonal and off-diagonal coefficients A Q( g 
depends on the intraband or interband interactions, re- 
spectively. In conventional materials with predominance 
of intraband scattering 3 F = A QQ Vp, so that the ef- 
fect of vertex corrections in Eq.(3) and (4) can be reab- 
sorbed in the definition of the transport scattering time 
r a = A aa /2T F , and the result (1) still holds, with renor- 
malized mobilities. Things are however deeply different 
in multiband systems with dominant interband interac- 
tions connecting e and h sheets, as in pnictides. In this 
case the largest elements in A a p are off-diagonal, leading 
to a mixing of the e and h characters and resulting in 
unconventional features, like a possible vanishing current 
3 F . In this situation, although the conductivities are still 
diagonal in the band index a (as one can show following 
the general derivation of Ref. [19]), Eq. (5) cannot be 
reduced to the Boltzmann-likc result (1). 

To investigate in details this issue, in the following we 
compute explicitly both Tp and 3 F in the representative 
case of pnictides, where the carriers in the h and e bands 
interact via spin- fluctuations (SF) exchange [1]. Accord- 
ing to neutron-scattering experiments [22], the SF spec- 
trum can be phenomenologically modeled with a stan- 
dard marginal-Fermi-liquid spectrum, 
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where \Q = Xo®/(T + 6) is the strength of the SF, 
oj s f = ujq{T + 0)/0 is their frequency scale and £t = 

£oV e /( T + ) is tnc AF correlation length, with 
Curie-Weiss temperature. Since the SF are peaked 
around the T — M nesting vector q = Q, the interac- 
tion mediated by such a collective mode will have a pre- 
dominant interband character. The crucial role of such 
interband retarded interaction has been already demon- 
strated for the understanding of several spectroscopic [23- 
25], thermodynamic[24, 26] and optical[27] anomalies of 
pnictides. However, to compute current vertex correc- 
tions the explicit momentum dependent of the bosonic 
spin spectrum (7), neglected so far in Refs. [23-27], must 
be taken into account. For the sake of simplicity we as- 
sume in the following that only interband scattering is 
present, neglecting any intraband coupling. The single- 
particle Green's function in each band is computed as 
usual by means of the Dyson equation G Q_1 (k, iuj n ) = 
iuj n — — S"(k, iuj n ), where the self-energy is given by 

S a (k,w n ) = 5 2 T]T X (q,^)G^k _ q, Wn - w,), (8) 

where uj n , uji are fermionic and bosonic Matsubara fre- 
quencies, respectively, and g is the coupling to the 
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FIG. 1: (a) T dependence of the renormalized currents (filled 
symbols), as compared to the bare velocities (empty symbols). 
We also show (dashed lines) the numerators of Eqs. (11)-(12), 
that fix the overall sign of the currents, (b) T dependence of 
the Hall coefficient Rh compared to the Boltzmann result (1) 
R H , computed with l/r Q = 2rJ. The units are fixed by the 
two-dimensional results divided by the interlayer distance d 
= 6.5 A. Dashed line: Rh obtained including also the effect 
of the FS shrinking, (c) Longitudinal resistivity as a function 
of T compared to p° = e 2 n a T a /m) 1 in the Boltzmann 



approximation with l/r a = 2Vp. 



bosonic mode x(q, iu>i) of Eq. (7). In Eq. (8) we ac- 
counted already for the nesting condition by considering 
only interband terms, so that the most relevant fluctua- 
tions are around q = 0. As far as the vertex corrections 
are concerned, following a standard derivation[18, 21], 
the current at ui = is computed as 

J«(k) = v"(k)+ 5 2 ]T [°° ^F(e)G^(k + q,e) 

xG^k + q, £ )Im X «(q, £ )J^(k + q). (9) 

where F(e) = coth(e/2T) - tanh(e/2T) and G R - A is the 
retarded/advanced Green's function. To compute the 
conductivities (3)-(4) we need the quasiparticle scatter- 
ing rates r^, and the dressed currents 3 F at the Fermi 
wavevectors k^.. By accounting for the most relevant 
contributions to the integrals (8)- (9) one can obtain 
an approximated semi-analytical expression for all these 
quantities [28]. The imaginary part of the self-energy is 
thus given by: 
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where w q = u> s f(l + £f-q 2 ). At the same time by intro- 
ducing the velocity and current projection along k, i.e. 



J a 



k and v a 



v F • k, so that v h < and v e > 0, 



J e 



(l-XheXeh)- 1 ^ +X he V e ) (11) 

{l-\ he X eh )- 1 {v e +\ eh v h ) (12) 



where the matrix A a p of Eq. (6) has been expressed in 
terms of the temperature-dependent coefficients 
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one finds [28] for the renormalized currents the expression 



, M ^ + (^ +q ) 2 l^ + q ' 

(13) 

By close inspection of Eqs. (10) and (13) we see that X a p 
increases as the nesting condition is approached, since the 
largest contribution to the integrals comes from vectors 
around q ~ kp — k p . Moreover, at high T where £t — 0, 
so that w q is independent on q, X a p vanishes and J a — 
v a . Conversely, at low T, £,t increases and only the value 
f5q = contributes to the integral (13), leading to X a/ 3 — > 
1 and large prefactors in Eqs. (11)-(12). 

To elucidate the effect on the transport of such a 
scattering mechanism connecting e- and h-like bands 
we will consider a set of parameters appropriate for 
electron-doped Ba(Fei_ a; Co a ;)2As2. In particular we take 
l/2m = 70 meV, and we choose = 90 meV and 

E m ax = 66 meV (i.e k F — 0.307r/a and kp = 0.377r/a) 
to reproduce the data at 7% doping, where long-range 
AF order is no more present and our model applies. 
For the SF we use [22] cj = 15 meV, 6 = 90 K, 
£o/a = 3.6 and g 2 XQ = 0-8 eV. Since the low-energy 
description (7) is not expected to hold any more around 
a scale of the order of the room temperature , we rescaled 
= <W /( T + 6)cxp(-T/T cuf ) to account for a fast 
decay of AF correlations above T cut ~ 300 K. Finally, to 
mimic the residual scattering by impurities at T = we 
added a constant (isotropic) scattering rate Tq = 4 meV. 
The resulting currents for each band as a function of tem- 
perature, as evaluated from Eqs. (10)-(13), are shown in 
Fig. la, along with the bare Fermi velocities. Two rel- 
evant features emerge. First, we find a strong tempera- 
ture dependence of both J e and J h , which deviate signif- 
icantly from their bare values with lowering temperature, 
due to the increasing scattering from SF. Second, in the 
e-doped case considered here, where \v e \ > \vh\, the dom- 
inance of XheVe with respect to Vh in Eq. (11) is reflected 
in a change of sign of J h at low T. These features have 
a striking effect on the Hall transport, shown in Fig. lb, 
along with the Boltzmann result (1) computed without 
vertex corrections. As expected, BP H is small and weakly 
T dependent, as due to the almost perfect cancellation 
of the contributions from the h- and e-like Fermi sheets. 
On the contrary Rjj has a strong temperature depen- 
dence and it can attain a large negative value at low T, 
where the e-like renormalized current |J e | 3> \ Jh\ domi- 
nates the transverse conductivity (4). At the same time, 
the effects of the vertex renormalization are less qualita- 
tively relevant on the longitudinal resistivity with respect 
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FIG. 2: (Color online) (a) Dependence of the vertex- 
correction coefficient A; le at T — 100 K on the Fermi- 
wavevectors in the two bands. Notice that Xhe decreases as 
one moves away from the nesting line kp — kp. (b) Rh as 
a function of the Co concentration x (which adds x electrons 
per Fe atom) The arrows indicate the corresponding values 
of the T = Boltzmann result (1), which has a negligible T 
dependence on this scale (see Fig. lb). On the right axes Rh 
is expressed in units of the inverse number of carriers per Fe 
atom, defined as 7ip e = 0.32 x 10~ 9 /\R H [m 3 /C}\. 

to the Boltzmann result, as shown in Fig. lc. Indeed, the 
dependence of <j% x in Eq. (3) on the sign of the renormal- 
ized current leads to a compensation between the vertex 
corrections in the h and e bands. Our results provide thus 
a consistent picture for both longitudinal and transverse 
transport, in good agreement with the experimental find- 
ings [3, 4]. For the sake of completeness we show in Fig. 
lb,c also the effect of the weakly temperature-dependent 
FS shrinking arising from the real part of the self-energy 
(8) [23, 27, 28]. As one can see, this is irrelevant on the 
Hall transport, while it contributes in part to the tem- 
perature dependence of the longitudinal conductivity, as 
discussed in Ref. [27]. 

It is also interesting to address the effects of dop- 
ing. Indeed, as we mentioned above and as we show 
in Fig. 2a, the absolute value of the A a( g coefficients de- 
creases as the Fermi wavevectors move away from the 
nesting condition kp = k F , realized at half- filling. We 
expect thus that the effect of the vertex corrections on 
the Hall transport will be less relevant by further in- 
creasing the Co concentration x. We investigate this 
issue by making a rigid-band shift of the chemical po- 
tential with doping, without changing for simplicity the 
microscopical parameters of the SF spectrum. The re- 
sulting Hall coefficient is reported in Fig. 2b, where 
we also mark with arrows the corresponding Boltzmann 
value R° H at each doping. As one can see, the low- 
T enhancement of Rh induced by SF decreases with 
increasing doping, and for x = 20% Rh almost coin- 
cides with R H ~ 1/ex, as found experimentally [3]. The 
trend shown in Fig. 2b, where SF spectrum is kept con- 
stant, is already in good agreement with the experiments. 



Nonetheless, one could also expect a decrease of the AF 
correlation length £p with doping, leading to a faster sup- 
pression of vertex corrections. This effect could explain 
the results in isovalent-substituted systems, as for exam- 
ple BaFe 2 (As,P/Ru) 2 [5, 6] or La(Fe,Ru)AsO[7], where 
the change of magnitude (or even of the sign) of the Hall 
coefficient should be attributed to a weakening of AF cor- 
relations, since no significant change on the FS pockets 
seems to occurs [29]. Finally, we notice that for a hole- 
doped system the overall temperature and doping depen- 
dence of Rh would be exactly the specular one: indeed, 
when the system is doped with holes, one has in general 
k F > kp, so that J h \ ;§> | J e \ at low T and the transverse 
conductivity will have a predominant hole-like character, 
in agreement with the experiments [8, 9]. Thus, the same 
mechanism accounts for the unusual Hall effect measured 
in pnictides both in electron and hole-doped compounds. 

In conclusion we analyzed the Hall effect in a multi- 
band model where carriers interact via the exchange of 
SF. We showed that when interactions have a predomi- 
nant interband character and connect carriers of opposite 
e/h nature, the currents renormalized by vertex correc- 
tions are dominated by the character of the majority car- 
riers. By evaluating this effect within a simplified two- 
band model and a phenomenological description of SF 
we were able to reproduce the main puzzling features ob- 
served experimentally in pnictides, namely a strong tem- 
perature dependence of Rh with a large absolute value 
at low T for weak (e or h) doping, and a more ordinary 
Boltzmann-like behavior at higher doping. We notice 
that the mechanism discussed here is quite general and 
robust: thus, an analysis based on a full self-consistent 
approach for the SF[18] within microscopic multiband 
models [10, 11] is expected to add only quantitative re- 
finements to the present results. An open question is 
instead the role of vertex corrections across the antifcrro- 
magnetic transition, where the Hall coefficient has been 
found experimentally to show an even larger T depen- 
dence. Even though this issue is beyond the scope of 
the present manuscript, it certainly deserves further in- 
vestigation to complete our theoretical understanding of 
transport properties in pnictides. 
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